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Abstract 

The braneworld model of Dvali, Gabadadze and Porrati (DGP) provides an intriguing 
modification of gravity at large distances and late times. By embedding a three-brane in 
an uncompactified extra dimension with separate Einstein-Hilbert terms for both brane and 
bulk, the DGP model allows for an accelerating universe at late times even in the absence of 
an explicit vacuum energy. We examine the evolution of cosmological perturbations on large 
scales in this theory. At late times, perturbations enter a DGP regime in which the effective 
value of Newton's constant increases as the background density diminishes. This leads to a 
suppression of the integrated Sachs-Wolfe effect, bringing DGP gravity into slightly better 
agreement with WMAP data than conventional ACDM. However, we find that this is not 
enough to compensate for the significantly worse fit to supernova data and the distance to 
the last-scattering surface in the pure DGP model. ACDM is, therefore, a better fit. 

1 Introduction 

It appears that the acceleration of the expansion of the Universe is now indubitable: it has 
been independently corroborated by measurements of type la supernovae 00 EI an d cosmic 
microwave background radiation observations by the WMAP satellite 0j- 

The simplest explanation for such an effect is the existence of a positive cosmologi- 
cal constant. Unfortunately, the estimates for its natural value are at least 55 orders of 
magnitude too large (see |U EJ for a review) . Another possibility is that the vacuum 
energy is zero, and the dark energy comes from the potential of a slowly-rolling scalar field 
|Hll3liniinilI3inilIllIlllISlinilIHl- Yet another alternative is to modify general relativity 
so that the effective Friedmann equation predicts cosmic acceleration even in the absence 
of dark energy. HH- A simple approach along these lines is to add an R n term to 
the Einstein-Hilbert action, with n < 0, which will affect the cosmological dynamics in the 
infrared. |2"T] 

The case we will discuss in this work is that proposed by Dvali, Gabadadze and Porrati 
("DGP") |22M23ll23] who suggest that the observed universe might be a brane embedded in a 
higher-dimensional space-time, with Standard-Model fields propagating strictly on the brane. 
Gravity propagates in the bulk, but on-brane radiative corrections to the graviton propagator 
result in there being induced an additional, four-dimensional Ricci scalar in the action. The 
cosmological solution to this theory with a five-dimensional bulk has been demonstrated to 
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possess a "self-accelerated" branch, as the universe approaches a de Sitter phase even when 
the vacuum energy vanishes in both the bulk and on the brane. |25j . 

By finding the equation governing the growth of radially symmetric perturbations on a 
cosmological background with zero cosmological constants (bulk and brane), Lue et al. [2£j 
conclude that the growth of structure is suppressed in a manner inconsistent with obser- 
vations and, in the linear regime, the DGP theory is a theory with a varying gravitational 
constant. 

On the other hand, Ishak et al. j22| have simulated data sets for CMB and weak lensing, 
assuming that the cosmology is pure DGP, and they have found that, if such data are 
analysed assuming that the cosmological acceleration is driven by a form of dark energy, 
rather than a modification of gravity, the forms of best-fit w obtained from the two data sets 
are inconsistent. 

Attacking the problem with more generality and a using a different methodology, we 
show that the quantity which is measured as Newton's constant for a small perturbation in 
the Robertson- Walker background is dependent on the energy-density content of the uni- 
verse and is, in general, not equal to the (truly constant) gravitational constant driving the 
Friedmann expansion. This effect results in a change in the evolution of gravitational poten- 
tials in cosmological models using linear perturbation theory: the matter power spectrum is 
somewhat altered and the integrated Sachs- Wolfe (ISW) effect is much weakened. 

It is well-known that the observed CMB anisotropies have less power on the largest 
scales than is predicted by the conventional ACDM model ^EHlEHj- As a consequence of 
the suppressed ISW effect, we find that DGP gravity provides a better fit to the temperature 
anisotropies observed by WMAP than does ACDM. However, DGP also predicts a somewhat 
more gradual onset of acceleration than that expected in ACDM, which is not as good at 
matching the observed Hubble diagram of Type la supernovae if we require that the distance 
to the last-scattering be fixed. Taken together, we find that although the low-multipole CMB 
data is significantly better fit by DGP, GR is preferred when both CMB and supernova data 
are considered for a flat universe. 



2 On-Brane Field Equations 

We start with a single 3-brane embedded in a five-dimensional bulk M . We will assume that 
all the standard-model fields are confined to the brane. Gravity propagates in the bulk and is 
a fully five-dimensional interaction, but, as first proposed in the model of Dvali, Gabadadze 
and Porrati j22], there is an on-brane radiative correction to the graviton propagator re- 
sulting in the induction of an on-brane, four-dimensional Ricci scalar in the effective action. 
Generalizing the DGP model slightly, we will allow for both a nonzero brane tension, A and 
bulk cosmological constant, A. No fields propagate in the bulk other than gravity. We thus 
write down the action: 



S 
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K + K~ 
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X + C 



SM 



(1) 



where g is the metric in the bulk, is the 5D Ricci scalar, ^R is the induced 4D Ricci 
scalar, \ parameterizes a vector field such that \ — coincides with the position of the brane 
everywhere (we assume that such a vector field exists); a particular choice for this vector 
field is n^, defined in <@J|. is the trace of the extrinsic curvature of either side of the 
brane; this term is the Hawking-Gibbons term necessary to reproduce the appropriate field 
equations in a space-time with a boundary. The appropriate energy scales are represented 
by n 2 — 8%M^ 3 for the true 5D quantum gravity scale and /x 2 = 8nM^ 2 for the induced- 
gravity energy scale on the brane. Finally, £sm is the standard-model lagrangian, with fields 
restricted to the brane. 
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The ratio of the two gravitational scales, the cross-over scale, 

^=2^ (2) 

was shown by Dvali et al. 22] to appear in the propagator of the graviton in the 5D- 
bulk DGP theory. It is at this distance scale that the potential owing to one-graviton 
exchange undergoes a transition from 4D to 5D behavior. Also, Deffayet [23 has shown 
that an FRW cosmology enters a self- accelerated phase when H ~ r^ 1 . If we wish to 
explain the cosmological acceleration using this phase, then we need to set r c ~ 1 Gpc, 
i.e. approximately the current Hubble radius. Dvali et al. (22] show that, at distances much 
smaller than the critical radius, when expanded around the Minkowski background, /x 2 enters 
the gravitational potential as the constant of proportionality, i.e. [j? — 8nG. This forces us 
to set M 4 = M P] = 10 28 eV and, therefore, M 5 = 100 MeV. 

Varying the action leads to the following equation of motion: 

( 5 )G^ = K 2 T^ , (3) 

where the indices run over all the five dimensions, from to 4. We write down the induced 
metric on the brane as 

q^v = SV - n^n v , (4) 

with n M a spacelike vector field with unit norm, normal to the brane at x — 0- We can then 
write down the energy-momentum tensor as 

Tpv = -Ag MV + S MV S{x) , (5) 
with the on-brane contribution described by 



2 (6) 

Here, is the Einstein tensor obtained by varying the usual Jd 4 x^R, restricted to 

exist on the brane, while is the energy-momentum tensor resulting from varying £sm with 
respect to q^ v . As already mentioned, the bulk is empty, save for a cosmological constant A. 

We then proceed by following the methodology first proposed by Shiromizu, Maeda and 
Sasaki [.304 : we carr Y out a 4 + 1 decomposition of the theory <PJ and calculate the effective 
on-brane field equations. This result was first obtained by Maeda et al. in (HI]. The detailed 
derivation is presented in Appendix El 

We use the Gauss equation to calculate the Riemann tensor on the brane: 

(4) iJ a /37i = ^R^ p<T q a ^ q"^ 5 + K a ^ s -K a 5 K 01 , (7) 

where K^ v :— q a ^q /3 iy V a nf3 is the extrinsic curvature of the brane. We can compute this 
curvature by assuming a Z2 symmetry of the codimension; using the Israel junction conditions 
for the jumps in the induced metric and the extrinsic curvature across the brane [22, we 
obtain 

2 <?^ = qtu -%v-- [%»] = (8) 

1 



2JT+ = [K^] = k I S M „ - -q^Sj . (9) 

We now contract appropriately and, after some manipulation, we obtain the equation for 
the 4D Einstein tensor: 



A A « 4 \ A ^\ 2 \ \-' 



Ak 4 

— + j^f^ - , (10) 
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where the tensors / M „ and E^ v are defined as 

f, u := ±AA^ - ~ AfA va + U^A* - l -A^j (11) 
£W = (5) C\ ±f3q %q^ , (12) 

with 

A^ (4) G,„^ 2 v- (13) 

-E 1 ^ is the projection onto the brane of the bulk Weyl tensor ^C a p jS (_L signifies an index 
contracted with f^ v describes brane terms that are quadratic in t^„ and G^ u \ hence- 
forth we will drop the superscript ^ and assume all quantities are four-dimensional unless 
otherwise specified. In Appendix El we explicitly define the cross terms appearing in / M „ as 
Ifiu and see equations ijfil)) — llfifil) for precise definitions. 

We can recover the standard Einstein equations from 11 Oil by sending the brane tension 
A to infinity, i.e. in the limit when the brane is 'stiff'. 

The difference between DGP gravity and the standard Randall-Sundrum result of (20] is 
the presence of the 4D Einstein tensor in Spu, resulting in additional terms containing ^G^ v 
explicitly: we now have an equation which is quadratic in the Einstein tensor, i.e. curvature 
can act as a source of curvature, leading to a potentially non-zero Ricci in the vacuum (as 
obtained, for example, by Gabadadze and Iglesias in j.i.'il l3"l] for their Schwarzschild-like 
solution in DGP). 



3 Cosmological Solution 

The quadratic nature of iflOjl makes this formulation computationally imposing. In addition, 
it was already noted by Shiromizu et al. in pHIl that the transverse-traceless component of 
E^ v contains the information about gravitational radiation coming off and onto the brane as 
well as the transition from 4D to 5D gravity, and its evolution is necessarily dependent on 
the state of the bulk: the equation of motion for does not close on the brane. 

Nevertheless, the large symmetry of a Robertson- Walker universe allows us to make 
progress. Choosing Gaussian normal co-ordinates, with the brane positioned at y = and a 
Robertson- Walker metric on the brane, our bulk metric becomes: 

ds 2 = -N 2 (y)dt 2 + A 2 {t, y) ll] dx i dx j + dy 2 , (14) 

and we are allowed to pick the normalization such that N(0) — 1 by rescaling the time 
variable; 7y is a 3D maximally-symmetric spatial metric. We also define the value of A on 
the brane: a(t) := A(t,0). In these coordinates, the 4D tensors are significantly simplified, 
with only zero entries in the columns and rows corresponding to the dimension perpendicular 
to the brane. We thus obtain for G^ v : 

G°o = -3(^ + 4) (15) 

*J ( 16 ) 
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with the dot representing differentiation with respect to the t coordinate. In addition, we 
assume that the brane is filled with a homogeneous distribution of a perfect fluid, such that 

= diag(-/3,p,p,j>,0) (18) 
p = wp . (19) 
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This allows us to obtain the 0-0 component of the quadratic tensor / M „ of (11 II : 



f° — 



1 ' " 2 



H p- 3 — + — 

or a z 



(20) 



The remaining issue is the tensor E^ v . Since we are dealing with an isotropic and ho- 
mogeneous universe, we must set it to be just a function of the time on the brane. Since it 
is traceless, it behaves just like radiation and decays as a -4 ; because of this property, this 
term is usually referred to as dark radiation (see, for example, the review of Maartens 
For the moment we will set: 

E° a = § , (21) 

although ultimately we will set this term to zero. 

Substituting the above into 11 Oil and replacing - with H, we obtain a quadratic equation 
for the energy density (or, equivalently, for H 2 ): 

K 4 2 /AK 4 K 4 / , k\\ 

i -r2 > 
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- 3 { 1 + W)[ H + ^) + T{ A + —) + -^°> < 22 > 

where k = — 1, 1, depending on whether the spatial hypersurface on the brane has negative, 
positive or no curvature. We can write this relation as a version of the Friedmann equation 
with modified dependence on the Hubble parameter, 



So, provided that the A and C terms remain negligible, and H 2 + 3> 2-^, the evolution of 
the scale factor in the early universe does not differ from that in standard FRW cosmology. 
As H decreases the evolution becomes non-standard; this can be demonstrated more clearly 
by solving {22j for H 2 : 



,M 4 At 2 k u 2 fi 4 p 2 , . . s k 2 . C 



H< = ^ + '- (p+ V-- 1+ 2^ r - i+ ^ + X)--A-- i , (24, 



with e = ±1 representing the two possible embeddings of the brane in the bulk (see [25])- 
The above result has been obtained previously obtained by Collins and Holdom [HEj and 
Shtanov |37| . 

We are going to follow Deffayet j2^] by naming the e = +1 branch as 'self-accelerated' 
and e = — 1 branch as 'non-accelerated'. The reason for this nomenclature becomes obvious 
in the case of flat bulk and a zero brane tension and no dark radiation, i.e. k, X,A, C = 0. 
Equation pHl then reduces to that originally obtained by Deffayet {ibid): 



,M 4 , , r, M 2 /m 4 , P 2 



H = 2 ^ + T p + 2 ^^ + T p - (25) 

As mentioned previously, this solution contains a cross-over scale above which the self- 
acceleration term dominates: 



K 2 



(26) 



2fi 

For e = +1, once /Lt 2 p/3 <C r~ /2, H approaches the nonzero constant r" 1 = 2/i 2 /« 2 , and 
the universe enters an accelerated de Sitter phase. If we wish to use this model to replace 
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the effects of the cosmological constant, we should set r c to be approximately the current 
Hubble radius. The non-accelerated branch with e = — 1 behaves like the usual Friedman 
universe, with H 2 tending to zero in the flat-universe case. 

Note that there are two distinct limits in which we can recover the ordinary Friedmann 
equation in the presence of a cosmological constant, 

H 2 = ^(p + Pvac)-4r' (27) 

where p is the energy density in everything other than the cosmological constant. One limit is 
to simply let k 2 — > oo, decoupling the extra dimension entirely and leaving us with p vac = A. 
The other is to set the brane tension to zero, A = 0, and take both k 2 and —A to infinity 
while keeping their ratio constant, yielding 



Pvac = y-^. (28) 

In our investigation of DGP cosmology, we will set the brane tension to zero while leaving the 
bulk cosmological constant A as a free parameter and calculating its likelihood as determined 
by the data. Then A = corresponds to "pure DGP," while A — > — oo corresponds to ordinary 
ACDM. 



4 Linearized Equations and the Potential 

In this section, we will derive the Poisson equation for a perturbation of the Robertson- 
Walker background and demonstrate that it deviates from the usual version: in its linear 
regime, the DGP theory on the brane is a varying-G theory. This derivation will assume 
that we are able to neglect all terms not linear in the gravitational potential. This range of 
validity of this assumption is discussed in 

We start off by introducing scalar perturbations to the metric. Since we are only going to 
be dealing with on-brane directions, we can use the 4D formalism in the conformal Newtonian 
gauge, parameterizing the perturbed metric by: 

ds 2 = -(l + 2<P(x,t))dt 2 +a(t)(l + 2${x,t)) llJ dx t dx 3 . (29) 

This allows us to calculate the Einstein tensor up to first order in the perturbations: 

G° = -3H 2 + 6H 2 <P - 6H& + ^V 2 <P + <D(<P 2 ,& 2 ) . (30) 

The above approximation is permitted provided <P, & <C 1. If we consider distance scales much 
smaller than Hubble scale, i.e. upon taking the Fourier transform, {k/a) 2 <P » H 2t P, H$>, then 
from the first-order terms we recover the potential for the flat Minkowski spacetime: 

(1) G° = !rV 2 <£. (31) 

a 1 

Alternatively, by including the Hubble flow, we can obtain the 0-0 component of the cosmo- 
logical evolution equation for the gravitational potentials in GR: 

1 ■ u 2 

— V 2 <P + 3H 2 V-3H&= ~P$ (32) 

where p is the average matter/radiation energy density of the universe, Sp is the deviation 
from this mean, and the fractional density excess is defined as 5 :— dp/ p. 
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In DGP gravity, the 0-0 component of the right-hand side of the modified Einstein equa- 
tion (1101 . first-order in <P or if', is: 

<" [RHS]» = -4„£ (A + p) + ^Hi - ^ (ff 3 * +fl'*) + ^({„ - 2»*)ff 2 

In calculating the above, we have assumed that the effect of the perturbations of the Weyl 
tensor, in this equation at sub-horizon scales are insignificant. Since we are dealing with 
a quasi-static situation, gravitational radiation is negligible. In addition, this tensor encodes 
the transition of gravity from 4D to 5D behavior. However, this effect only occurs at a 
length scale determined by r c , which, as shown in the data fits of ii6.2l is much larger than 
the horizon scale even today, let alone in the past. As such, it is unlikely to have any 
significant effect on the quantities under consideration. 

Performing some algebraic manipulation and substituting for H 2 from jSH) we obtain the 
DGP equivalent of ®: 



1 a 1 
-^V 2 <£ + 3H 2 & - 3H$ = -—p 
a 1 2 



(34) 



If we apply similar approximations to those that led to l(3ll . i.e. assuming that <C 1, 
drop terms containing H by considering scales over which the Hubble flow is negligible, we 
obtain: 
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(35) 



This form of lj35|l signifies that the weak-field approximation in a Robertson- Walker back- 
ground is a varying-G theory, with the effective Newton's constant dependent on the average 
energy density of the Universe: 
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(36) 



Note that, in the numerical solutions described below, we do not neglect terms containing 
H; we are taking that limit here purely for expository purposes. 

The large-scale evolution of the universe is always driven by p3|). i.e. the energy scale 
encoded in /i 2 . This is the parameter which drives the conditions during, for example, 
nucleosynthesis. However, at least in part, the evolution of structure is driven by the effective 
Newton's constant, i.e. equation l.ifill : this is generally true at late times for diffuse clouds 
of gas. This insight will provide a constraint for some of the parameters of the theory. 

As we can observe from equation i|36|l and assuming that we are in the self-accelerated 
branch of the solution (e = +1), the value of the effective G e ff varies from one to, at most, 
two times the underlying Newton's constant, with the increase occurring in the late universe 
as p — > 0. (A related effect has been observed in the presence of Lorentz- violating vector 
fields [38 , which result in a cosmology where Newton's constant differs from the constant 
relating energy density to the Hubble parameter in the Friedmann equation.) In the limit 
described at the end of sectional in which we take A — > — oo, the time-dependent piece of 
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G e ff goes away, and we obtain G e g — /i 2 /87t = constant, so that ordinary ACDM is indeed 
recovered. 

By considering the i-j equations of motion in GR, we can get the an equation relating 
the two potentials <P and to the anisotropic stress of the cosmological fluid, ir: 

^(£ + ^)=/Z 2 7T (37) 

Since the anisotropic stress is negligible in models with no neutrinos, we can set <P = 
The analogous calculation in DGP is a little more complex; however, it yields similar results: 
without a source of significant anisotropic stress in the cosmological fluid, we can set <P = —fy . 



( K 2 ( 1 , , 3a \\ k 2 , 

k 2 ( 1 _. , 3a \ 2 

= 6^( v A+ 2^ (1 - 3u;) -^J^ 7r - 7r£ 

tte is a new term and is the anisotropic stress encoded in the the Weyl tensor E^ v , i.e. a 
result of off-brane effects, such as graviton evaporation into the bulk and any gravitational 
waves going off or coming onto the brane. We will set this to zero in the subsequent analysis. 



5 Range of Validity of Linear Regime 

In section^] we assumed that we can use the linear approximation to the modified Einstein's 
equation iflfljl . It is important to verify this explicitly, since the large magnitude of the 
coefficients of the quadratic terms in ifTOj) . might lead to their dominating over the linear 
terms. 

After a systematic review of all the quadratic terms in the expansion of IKH , we conclude 
that the quadratic term which is most likely to be large arises from the (G^) 2 term and is 
of the form 

£ (V 2 ^) 2 (39) 

Assuming that the Poisson equation is approximately valid, despite the evolving cosmological 
background, we have V 2 ^ ~ A* 2 p<5, and using the definition of r c , ll21it . leads to the condition 
that for the linearization of DGP to be valid 

r 2 cf i 2 p5 < 1 (40) 

If we assume that the density perturbation is in the form of a spherical top hat with a 
characteristic size D, we can write down the mass of this object as M ~ D 3 pS. Finally, the 
Schwarzschild radius in four dimensions is rs ~ ^ 2 M, giving us: 

r 2 c r s < D 3 (41) 

We have thus recovered the result first discovered by Dvali et al [HHI, that there is a new 
scale in DGP theory, 

U = (r c 2 r s ) 1/3 . (42) 

At small distances (between the Schwarzschild radius and r*), gravity behaves essentially as 
in 4D GR, because the quadratic terms dominate the modified Einstein equation iflOj) . If we 
set A = X = 0, in this regime we are looking for the solution to 



(43) 
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with the obvious solution being G ll v = Ixt^ , i.e. the usual 4D Einstein equation. On inter- 
mediate scales (between r* and r c ), gravity is described by a scalar-tensor theory, consistent 
with our description above in terms of a time-dependent gravitational constant. It is in this 
regime that the linearized DGP equations are valid. 

As we will see below, as a perturbation of fixed comoving size expands along with the 
universe, it typically goes from being less than r* to being greater than r*. Each mode is 
therefore described by 4D GR at early times, and later on by linearized DGP. 



6 Cosmological Simulations 
6.1 The Simulation 

In order to test the effects on cosmological observations of modifying gravity, we built a 
simple cosmological simulation containing only dark matter and radiation and modeled their 
evolution according to linearized equations - linearized DGP or ordinary 4D GR, depending 
on the regime a given mode is in. The model produces as outputs the matter power spec- 
trum and the contribution of the Sachs- Wolfe effect (both integrated and non-integrated) 
to the radiation power spectrum. The difference between the matter power spectra is slight; 
however, the impact on the ISW is significant, resulting in a large reduction of power at low 
multipoles. 

We compare the results of the concordance ACDM model to that of DGP cosmology. 
The concordance model is defined by Q m = 0.27, O m /S2 T — 3234 and Q\ + £2 m + Q r = 1. 
For the DGP cosmologies we use the same matter to radiation energy ratio, and choose the 
self-accelerated embedding, e = +1. We keep fi r h 2 constant. As our aim is to explain the 
acceleration of the expansion of the universe using effects arising from modified gravity, we 
have set the brane tension A (which is equivalent to the cosmological constant in ACDM) to 
zero in all our considerations. For simplicity we have also set C and k to zero, i.e. we are 
assuming a flat cosmology. We keep the bulk cosmological constant A as a free parameter. 

We parameterize the DGP models by rewriting ll24T) as: 



H 1 1/1 

H 2 = 2/32 + Qma Z + Qra i + p V 4^2 + ^ma~ 3 + ^i-a- 4 - Q A (44) 
where H$ is the Hubble parameter today, 

(3 := r c Ho cross-over radius in Hubble units 



J? m + J? r := (^po/SHq = 1-/3 1 \J 1 — S7a energy density 

fl A k 2 A/6Hq dimensionless bulk cosmological constant 

S7 r = J? m /3234 fractional radiation density 

f2 m is the contribution of matter to the total energy density of the universe (the remainder 
coming from the DGP curvature). Realistic DGP models will have f2 m ~ 0.3, just as in 
ACDM. The energy scale for the bulk cosmological constant is A « (10 -8 eV) 5 .!?^. 

A significant issue in modeling DGP is deciding on where the the transition between the 
Einstein and the linear DGP regimes occurs and how to implement it. The simulation was 
built to switch from GR to linearized DGP instantaneously at the point in evolution defined 
by the scale r*. We will apply DGP gravity when the scale of the perturbation is 

D 3 ~rl = 2GMr 2 c = ±^D 3 r 2 p8 (45) 

The perturbation size is eliminated from the relationship, yielding for the fractional density 
perturbation at which GR transitions to linear DGP: 

24 



p 2 r 2 p 



a 3(l+w(a)) ( 4 g) 
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where w(a) is the effective equation of state parameter for the fluid comprising the universe 
and a is the scale factor; S is dependent on the initial conditions, and thus is a function of 
the mode under consideration. Thus, for matter and radiation, the transition occurs from 
GR to DGP and occurs only once in the evolution of the perturbations. This transition has 
been plotted on Figure We can determine the success of the splicing between GR and 
the linear DGP regimes by comparing the value of the effective Newton's constant at the 
transition point. We found that the transition occurs at around z = 20 for modes of the size 
of the horizon today, with G e g « 1.02G, with this transition occurring at larger redshifts for 
higher modes and G e ff even closer to 1. 




Figure 1: Value of scale factor at which the gravity driving evolution transitions from GR 
to linear DGP. Solid line is for Minkowski bulk = 0), dashed for Qa = — 1. Shape of 
transition line is a processed dark-matter power spectrum. 

In Figure we plot the evolution of the Newtonian gravitational potential <P for different 
modes in both ACDM and DGP. Each mode diminishes in amplitude when it first comes 
into the Hubble radius (more so during radiation domination), before settling to a constant 
amplitude. In ACDM, linear modes once again decrease in amplitude when vacuum energy 
begins to dominate over matter. In DGP, in contrast, the effective gravitational constant 
increases just when the universe begins to accelerate, leading to an increase in the amplitude 
of $ that can be appreciable for certain wavelengths. This behavior is also reflected in the 
transfer function for dark matter perturbations, plotted in Figure which shows a slight 
enhancement in DGP over ACDM. 

6.2 Constraints from Expansion History 

Before considering details of CMB anisotropies, we turn to two sources of constraints on 
the expansion history of the universe: the Hubble diagram of Type la supernovae, and the 
distance to last scattering as measured by WMAP. These imply a tight relationship between 
the two free parameters in the DGP model: /3 = t c Hq and Qa- We will then calculate CMB 
anisotropies in models that obey this relationship. 

For the supernova data, we used the Riess et al. Gold SNe la data set @H] (156 super- 
novae) and searched for the parameters minimizing x 2 . In the DGP scenario with both (3 and 
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log(a) 

Figure 2: Evolution of Newtonian potential in GR and DGP (GR — solid lines, DGP — 
dashed). In DGP, the growth of the effective Newton's constant leads to a wavelength- 
dependent growth in the potentials at late times, as opposed to the decay observed in ACDM. 
(ACDM: concordance model; DGP: Q A = 0,/3= 1.38) 
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Figure 3: Transfer functions in ACDM and DGP (f2 A = 0, (3 = 1.38). The transfer function 
is defined here as ratio of size of initial perturbation in <P to its final value, renormalized to 
1 at large scales. In GR, the growth rate for the potential during dark-energy domination is 
independent of k, so the quantity shown here is equivalent to the usual transfer function. For 
DGP, the growth rate depends on k. 
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Qa free, the optimization routine chooses very large and negative values for Qa, correspond- 
ing to the ACDM limit as discussed at the end of Section 03 The supernova data, in other 
words, prefer a very large and negative bulk cosmological constant, reducing the observable 
physics to GR. This is the result of the fact that ACDM fits the supernova data slightly 
better (see Table 0. In terms of the effective equation-of-state parameter w, the supernova 
data prefer w w — 1 or even a little less, while pure DGP corresponds to w e g ~ —0.7 today. 

If we fix Qa = 0, supernova data fit best when (3 = 1.20, implying Q m — 0.20. The 
details of the fit are presented in Table 

The WMAP experiment has determined, to a high level of precision, the distance to the 
last-scattering surface as d\ s — 14. 0^3 Gpc = 3.32+°;^ -^o" 1 E9- This restriction needs 
to be included in the likelihood calculation. Whereas for ACDM the two sets of data are 
consistent, they are not so for pure DGP. A good fit for the CMB distance in DGP with 
Qa = requires a higher (3 = 1.66, implying Q m = 0.40. 

This is very different from the requirements of the supernovae. Putting the two sets of 
data together results in the conclusion that owing to the inconsistency between parameters 
preferred by the SN and CMB distance data, the overall fit to pure DGP assuming a flat 
cosmology {Qa = 0) is much worse than for concordance ACDM. 



Scenario (3 = r c Ho Qm X 2 P er d.o.f. Confidence 



SN DGP 


1 26 +u - Ui 


0.20 ±0.01 


1.15 


9% 


SN ACDM 




0.30 


1.14 


10% 


CMB dist. DGP 


1 66+ - 08 

1 - DD -0.02 


40+ 002 






CMB dist. ACDM 




0.29 






Total DGP 




0.28 ±0.01 


1.21 


4% 


Total ACDM 




0.29 


1.14 


11% 



Table 1: A set of best-fit parameters for yICDM and DGP cosmologies with Q A = 0. 'SN' 
represents fits to just SN data; 'CMB' are fits to the distance to the last-scattering surface; 
'total' combine the two data sets. For ylCDM, the two fits are consistent, for pure DGP they 
are not, resulting in a significantly worse overall fit. 

Using both the SN and WMAP data, with both Qa and (3 free, the maximum likelihood 
for the DGP cosmology is attained for Qa — —7.3 and (3 — 4.1. However, since the range of 
admissible values of (3 increases for more negative Qa, the likelihood for just Qa (marginalized 
over (3) increases monotonically as Qa attains lower values. Figure^lshows the likelihood for 
(3 given a particular value of Qa , while figure shows the likelihood for Qa marginalized 
over the supernova absolute magnitude and (3. From the expansion history alone, ordinary 
GR (corresponding to Qa — > — oo) is preferred. 

6.3 Simulation Results 

The matter power spectrum for DGP is slightly different than that for ACDM. We have found 
that there is excess power at large scales and a deficiency of power at low scales. Results 
are shown in Figures and 03 This change is a result of the different late-time evolution 
of the gravitational potential and the change in the rate of growth of density perturbations 
associated with it. 

Using the sets of parameters presented in Figure 21 we computed the ISW effect for 
the DGP model. We have found it to be significantly reduced at low multipoles and to 
have a (very) slight excess in the power for £ above 20 as compared to the concordance 
ACDM model. Making the bulk cosmological more negative restores the behavior observed 
in ACDM. A positive Qa significantly increases the ISW effect. The results are shown in 
Figured 
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Figure 4: Plot showing the one sigma (dark color) and three sigma (light color) range for 
best-fit values of (3 for given values of Qa for CMB distance (green), SN (blue) and combined 
(red) For Qa close to the preferred values for the two data sets are significantly different, 
leading to a poor overall fit. As Qa — * — oo the preferred parameter spaces increasingly 
overlap. In this regime DGP is indistinguishable from GR. The table presents values of best- 
fit /3's for a selection of Qa and the x 2 ' s of the respective fits to combined CMB-distance 
and SN data. It can be clearly seen that a positive bulk cosmological constant is strongly 
excluded. 



6.4 Impact on CMB 

From the fits performed in M6.21 we obtain a range of values of (3 for each Qa which satisfy 
the constraints of the supernova data and the distance to the last-scattering surface. In 
our subsequent analysis and simulations we pick as (3 the central values of the likelihood 
distributions for each Qa • A selection of Qa and (3 pairs is presented in the table contained 
in Figure Q] 

By requiring that the distance to the last scattering surface is effectively fixed and by 
assuming the same radiation-density-to-matter-density ratio as in the concordance model, 
we ensure that the part of the CMB spectrum resulting from plasma oscillations remains 
unchanged, despite the different theory of gravity. We are thus able to take the ACDM 
output of CMBFAST and 'replace' the ISW part of the power spectrum with the new DGP 
calculation, provided we properly take into account the cross-correlations between the ISW 
effect and the other contributions to the CMB power spectrum. 
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Figure 5: Comparison of dark-matter power spectra for ACDM (solid line) and DGP (Qa = 
0,(3 = 1.38, dashed). The spectra have been normalized to unity at k = 10~ 3 /i Mpc -1 . There 
is excess power at large scales and a power deficiency at low scales. 
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Figure 6: Ratio of dark-matter power spectra for DGP (Qa = 0,/3 = 1.38) and ACDM. 
Both spectra are initially normalized to unity at k = 10~ 3 h Mpc -1 . 
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Figure 7: Contribution to Ct from ISW effect as calculated for ACDM and a range of DGP 
cosmologies. The effect is significantly reduced at low multipoles in DGP. Making the bulk 
cosmological constant more negative brings the effect back to ACDM levels. Making Qa more 
negative brings the DGP power spectrum towards the ACDM one. 

Since the low multipoles are dominated by the SW and ISW effects, we can compute 
both the ISW and SW effect in all cases and then assume that the correlation between the 
two is equal to the correlation between the ISW contribution and the rest of the signal. The 
procedure is explained in detail in Appendix [Bj 

The results of performing the above procedure for models for DGP cosmologies with a 
range of S7a (with (3 as implied by SNe la data and d\ s as from WMAP) are presented in 
FigureEl The total signal strength at low multipoles is significantly reduced. The likelihood 
for Qa as implied by the fit to the low-multipole WMAP data alone is presented in FigureEl 
For this data alone the maximum likelihood is found at Qa — +0.06, with a bulk cosmological 
constant higher than that strongly excluded. Thus, the CMB data alone slightly prefer pure 
DGP to ACDM. 

Finally, we combine the likelihoods obtained from the SNe la data, and the CMB data, 
presented in Figure El to obtain an overall likelihood distribution for Qa- This has been 
presented in FigureEB Taking all experimental data together does not change the conclusion 
that ylCDM is preferred to DGP: the better fit to low-^ multipole data of WMAP is not 
enough to compensate for the inconsistency between fits to CMB distance and supernova 
data. 

7 Conclusions 

We have performed a projection of the equations of the Dvali-Gabadadze-Porrati modified 
theory of gravity in 5D bulk (Minkowski, de Sitter and anti-de Sitter) onto a 4D brane 
embedded in it. We have rederived a cosmological solution to the theory and have derived 
equations governing the evolution of linear perturbations. We find that the theory governing 
the linear perturbations, in certain regimes, is one where G is not a constant, but is dependent 
on the average energy density in the universe. 

Using the new equations, we built a simple cosmological simulation containing radiation 
and dark matter, driven by DGP gravity. We have discovered that, in DGP cosmologies, the 
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Figure 8: Output of CMBFAST for concordance ACDM (blue solid line) and its modification 
for pure DGP {CI a = 0, red dashed line). Shaded are represents cosmic variance for pure 
DGP. WMAP measurements have been plotted in black. DGP reduces power at the lowest 
multipoles, bringing the power spectrum into slightly better agreement with measurement. 



0.07 



0.06 - 



0.05 - 



g 0.04 - 
'a; 

^ 0.03 - 



0.02 - 



0.01 - 



i 1 1 1 1 1 

-20 -15 -10 -5 5 

n A 



Figure 9: Likelihood for a range of values of CIa as implied solely by data from low-multipole 
WMAP measurements. The preferred value is Qa — +0.06 with the upper half-maximum 
lying at lying at +0.3. Normalization is arbitrary. 
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n A 

Figure 10: Likelihood for a range of values of Qa as implied by data: dotted line includes 
only CMB distance and SNe la. Solid line is the modification to the likelihood caused by 
adding constraints from low-multipole WMAP data. Despite the fact that WMAP data 
prefers pure DGP, the preferred choice for the value of the bulk cosmological constant is high 
enough to be excluded by the other data. ACDM is clearly preferred. 

Newtonian potential <P exhibits a period of growth at late times, prior to its decay; this is 
in contrast with GR, where the potential decays once dark energy dominates. The precise 
details of this effect are a function of the wavelength of the perturbation, and lead to an 
altered transfer function and a changed dark-matter power spectrum, with slightly higher 
power at large scales. 

The impact of the change in the evolution of the potentials can be seen through the 
Integrated Sachs- Wolfe effect. We simulate the CMB anisotropy and demonstrate that in 
DGP cosmologies the ISW effect is significantly weaker at low multipoles. 

We constrain the parameters of the theory through a fitting procedure using supernova 
data and the WMAP results and perform a calculation of the likelihood function for the pa- 
rameter space. We find that assuming a flat cosmology, pure DGP with no bulk cosmological 
constant is significantly worse when simultaneously fit to supernova data and the distance 
to the last-scattering surface. DGP fits better to low-£ multipole data from WMAP, owing 
to the reduced ISW effect. However, this does not compensate for the poor fit to the former, 
leading to the conclusion that GR and /1CDM are preferred by all the data available taken 
together. 
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A Derivation of On-Brane Field Equations 

In this appendix we will replicate the work of |3()U31j . deriving the projected, on-brane field 
equations for DGP gravity in 5D bulk. 

We start off with Einstein's equation in a 5D manifold M 

G pv = n 2 T^ . (47) 

The Greek indices run over all the dimensions, i.e. from to 4. 

We want to put in a 3-brane, B, (we will denote its position by invoking a spacelike vector 
field x m the neighborhood of the brane, such that \ — will coincide with the position of 
the brane. We want then to find the effective equation of motion for gravity on B itself. 

Let us first take care of the LHS of l(T7j) . We will first define the induced metric on the 
3-brane, q^: 

%lv = g pv - n^i u , (48) 
where g pv is the metric on M, and is a spacelike vector field in M, with unit norm, 
which is normal to the brane at x = 0. Now we invoke the Gauss equation to calculate the 
Riemann tensor on the 3-brane 

^R a W = (5) #%a<Z>> P 7 <f 5 + K%K PS - K a s K 01 , (49) 

where K pu :— q a fl q^ 3 iy W a np is the extrinsic curvature. 

In order to simplify our notation, we are going to define the perpendicular subscript to 
signify A pl _ := A pv n v . 

Contracting a and 7 in ijiflf we obtain the 4D Ricci tensor: 

w Rps= {S) R v<7 q v 0%- WR^-^tfpfs + KKps- K a K Sa . (50) 
and, then, the 4D Ricci scalar: 

&R:= ^R ltu g^= {5) R^q^- ^B^^qT + K 2 - K»K\ (51) 
= &R ±± - ^R ± p /+ ^R^^+K 2 -K a[i K a P . 

=0 

An observer on the brane will only be able to measure the 4D Ricci tensor and scalar, 
and not their five-dimensional counterparts, so will naturally define the 4D Einstein tensor 



WGps := {i) R 0S -\qp S (4) i?= (52) 

= < 6 >jwy*- (B) «w*-^^( (5)i2 - (5)R ^- {5)R \± a ) 

+ KK fjS - K p a K Sa - X - qs)8 (K 2 - K al3 K al3 ) . 
We can simplify this by: 

^R ±1 _ = g a(i rf WRp^x = 0°V {b) R\ ±a = (5) i?%/ ■ (53) 

and also expanding the Riemann tensor in terms of the Riccis and the Weyl tensor in 5D 
gU (3.2.28)]: 

q a y u n p n° (5) iW = q a ^ v n p n° Q (g pa &Rp a - g pP ^R aa - g aa ^R Pp + (54) 

+g a p {5) Ra P ^j - (g P agp a - g P pg<7 a ) {5) R + (5) C paal3 ^j = 

= I ( {5) R a pq°y„ + q pv (5) Hxx) - ^q pv {5) R + ^C^q*^ . 



= -E u 
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We can also manipulate the 5D Einstein equation I|47J1 to obtain: 



^R=-\n 2 T (55) 

{5) R P a = k 2 (r pa - \g P aT^j (56) 
R±x = k 2 (t ±± - hf\ . (57) 

Using all of the above in il 5 2 fl . we obtain: 

(4) G,v = (5) G pff q/g/ + {5) R^q^ + KK^ - K/K vp - l -q pv (K 2 - K al3 K af3 ) - 

- I (5) R aP q a ^v ~ (5) R±± + {5) R- e pv = 

= ^- (T P aq^q\ + (t ±± ~ \ T )) + KK ^ ~ K/K vp - 

~ \%v {K 2 - K^K aP ) - . (58) 

Now let's take care of the energy-momentum tensor. It has contributions from the bulk 
(which we will limit to just a bulk cosmological constant, A) as well as the brane contributions. 
One of the contributions in the brane is the radiatively induced DGP 4-gravity term; the 
action here is just the usual ^~ 2 S(x) J d 4 £ ^ R, restricted to the brane, the variation of which 
gives the 4D Einstein tensor. The stress-energy tensor has the form: 



T 



-A 9pLV + (-\q p „ + r M „ - ii- 2 ( 4 >G^) S( X ) • (59) 



We are going to perform our calculation just off the brane. The energy-momentum 
tensor there consists of just the bulk cosmological constant. The extrinsic curvature can be 
determined from the Israel junction conditions, one of the formulations of which allows us 
to compute the jump in the metric and extrinsic curvature across a thin shell of non-zero 
momentum-energy |32j : 

Q-tu - %v =■ \%A = (60) 
[K^] = -n 2 [s pv - p^s) (61) 

Assuming that the universe is symmetric about the brane allows us to obtain the values of 
the extrinsic curvature explicitly, Kf v = —K~ = \[K^] and K + = ^n 2 S. Substituting 
into (l58l) and making 4D subscripts implicit: 

G pv = ——Aq^ v — — S (s - -qfj, v Sj — (s^* - ^QffSj \ S va — -q va 

We then substitute for S^„ from its definition and, after a bit of algebra, obtain: 
I 6/? J = ~ V~2~ + J ^ + + K + -j7*» + -pfa - E ^ ■ (63) 
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where the new tensors are quadratic in G M „ and t^ v and are defined as: 

7r MV = - ~T„°T„ Q + ~<? M „ (V^T^ - ^T 2 ^ (64) 

7m w = To^'^'W ~ ,^G va + —9^11/ ( G a pG a — ) (65) 



£nv = ~jx (tG^ + t^G) + - (r^ a G ai/ + G^t^) - -q^ ( T al3 G at3 - -tG J (66) 



12 

^ - C^pif^ . (67) 
These quadratic tensors are related to / M „ of 1111) by 

fixu = ^K/iu + fJ 2 ^ + JfJ,U • (68) 

B Modifying CMBFAST data 

This appendix explains how CMBFAST data was modified to reflect the expected DGP 
signal. 

Our cosmological simulation was used to calculate the contribution to the CMB power 
spectrum from the ISW and SW effects, as well as the total of the two (taking into account 
the cross-correlation) . We then defined the IS W-SW cross-correlation factor as 

/-"ISW-SW /-(total pSW-SW r<ISW-ISW 

R := , G < = = ^ ~ C * = (69) 

r'isw-iswr'sw-sw o /nsw-iswr^sw-sw 

The R thus defined was then assumed to be identical to the cross-correlation factor 
between the ISW effect and the total CMBFAST signal. This is a very good approximation 
at the lowest multipoles. With this assumption we can compute the non-ISW part of the 
power spectrum by solving for C^ ISW in 



^CMBFAST ^~ISW + 2 / Q~ISW QlSVf-ISW R _|_ ^ISW-ISW ,jq\ 

Finally, to compute the power spectrum for DGP cosmologies, we take the values of R and 
£<isw-isw ca i cma ted for the particular model and use them in equation ijTfljl . This output is 
presented in Figure 03 
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